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Abstract. We discuss a Brans-Dicke model with a cosmological constant, negative value of 
the w parameter and an arbitrary (in general non-vanishing) scale factor at the Big Bang. 
The Friedman equations for a flat universe are considered. The current observational values 
for Hubble constant Hq and deceleration parameter qo play the role of initial conditions. 
We follow the approach of [1] in order to solve field equations analytically. In Ref. [1] only 
positive values of w were considered, we extend the study to a complete set of possible w 
values. Our main result is that the scale factor (during it's evolution back in time direction) 
may not vanish, unlike in the standard ACDM case. In other words, the considered model 
demonstrates a cosmological bounce instead of the initial singularity. The famous formula 
(4.6), that leads to the bounce, is valid only for the dust-filled universe with p = and, 
therefore, is not adequate for the Early Universe hot stage when the bounce happens. So, 
our results are qualitative in nature and must be used to obtain initial values for the hot 
stage of the Universe. 
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1 Introduction 

The cosmic acceleration is now a well observationally established fact [2-7], however it's 
physical reasons remain open. So, the construction of an adequate cosmological model with 
a minimal number of ad-hoc hypotheses is a hot topic in modern cosmology. The simplest 
and hence the most studied model for the moment is cold dark matter one with cosmological 
constant {AC DM). Providing a good quantitative agreement with observational data, this 
model however does not explain the nature of the dark matter and dark energy. Another 
weakness of AC DM is the absence of explanation of the smallness of A value if it is assumed 
to be the so-called "vacuum energy". All these arguments lead to the idea of a dynamical 
theory of dark energy creation (see for example [8] ) . the most widely discussed candidates are 
quintessence (a slowly rolling scalar field [9]) and higher order curvature gravity (including 
so-called f(R) gravity models [10]). 

Brans-Dicke model (BD) is one of the first gravity models with a scalar field [11]. 
It was suggested in 1961 and contains an additional parameter uj whose value has to be 
determined by observational data. Large values of w mean an important contribution from 
the tensor part (Ricci scalar) , smaller values of w mean an increasing role for the scalar field 
contribution. In the limit w — > oo BD theory leads to General Relativity (GR). In BD model 
the value of Newton's constant is proportional to the inverse scalar field (G ~ l/</>), proving 
additional coupling between the model parameters. The most accurate limit on w comes 
from Cassini-Huygens mission data on post-Newtonian parameter 7 and is \ W\ > 50000 [12]. 

BD theory is the most natural GR extension. It is interesting because, firstly, this 
model could be the low energy effective limit of grand unification (and super unification) 
approaches (from latest LHC data (see [14]) this possibility is not completely closed yet). 
Second, because the scalar field in BD theory can be reinterpreteted as a dilation fiels in 
string theory. Finally, because BD model is the simplest GR extension and is useful to 
investigate any supertheory, so as to gauge the difference with GR [13]. In addition, BD 
gravity is widely used in cosmology as one needs a scalar field for inflation and such a filed 
is anyway necessary in the BD model. A large set of inflationary models [15-17] is based 
on BD gravity and more generic scalar-tensor approaches. Brans-Dicke theory is also closely 
related to the widely discussed f(R) gravity (see for example [18]). 
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It is necessary to underline that there is no accelerated expansion in the standard version 
of the BD model, so, one has to study it's extended versions. One of the most common 
extensions is a scalar field potential addition. As the accurate shape of this potential is not 
known yet [20], one can consider a A-term as the effective contribution instead of a potential 
(so we obtain BDA model). The explanation of the A-term smallness in the framework of 
BDA is possible and was suggested in [21]. In Ref. [22], with the help of scalar field in BDA 
context, a dark matter hallo around galaxies is modeled. 

An analytical accurate solution for Friedman equations in the BDA model was also 
obtained in Ref. [1] where positive values of u and initial conditions for the scale factor in 
the form a{t m i n ) = 0, where t m i n is Big Bang time, were considered. Partial solutions in this 
model with scalar field power dependence versus the scale factor were presented in [23-26]. 
Vacuum solutions were obtained in [27-29]. Some papers discussed a A-term dependence as 
a function of the scalar field (for example, [30]). Numerical integration and stability analysis 
of .BDA+matter solutions were carried out in [31]. Big Rip solution in BDA is discussed 
in [32]. The analytical solution in the pure BD model with negative w, avoiding the initial 
singularity, was obtained in 1973 by Gurevich et al. [19]. There is no cosmic acceleration in 
the Gurevich et al. solution, so, nowadays one has to extend BD theory to include cosmic 
acceleration. In this paper we therefore study the Einstein-Friednman equation solutions 
in BDA theory for w < with a scale factor with initial value a(t m i n ) = a m i n . Generally 
o-min 7^ 0. Friedman equations are studied for a flat Universe. We take the current values for 
Hubble parameter and its derivative (deceleration parameter) as their initial conditions. In 
our study, we use the approach suggested in [1]. We would like to point out that only positive 
values of w are considered in [1], so, our solution with w < represents a new branch. As 
opposed to the standard ACDM-model, in the considered case, the scale factor never vanish 
during backward time evolution. A so-called "bounce" (a snap back from the minimal value 
of a m i n ) corresponds to this situation. The expression 4.6 leading to the bounce is obtained 
for a cold Universe with p = and is not valid for a hot Universe. Therefore all the values 
in the bounce region are only qualitative estimates for the initial values for the transition to 
the hot stage. 

This paper is organized as follows: 
in section 2 we discuss the choice of the space-time metric and the corresponding field equa- 
tions; in section 3 the initial values for cosmological parameters are obtained; in section 4 
we obtain an analytical solution with a bounce for a dust-filled universe (p = 0); section 5 
contains a preliminary discussion of the results of section 4; in section 6 we explore the case 
of a ultrarelativistic state of matter (hot phase); section 7 is devoted to the conclusions. 

2 Field equations 

TThe Friedmann- Robertson- Walker (FRW) metrics reads 1 : 



ds 2 = dt 2 - a(t) 



2 f dr 2 
1 — kr 2 



+ r 2 dVi\ 




where k = 0, ±1. 

The action of the BD-A theory can be written as: 



S = — d 4 x^g $(R + 2A) - -g^d^d^ + 16vrL, 




matter 



(2.2) 



Here and below we set light speed c = 1. 
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Here w is the BD theory parameter, $(i) is the scalar field, A is the cosmological constant 2 . 

Variation of the action with respect to the metric g^ v and the scalar field $ gives the 
following field equations: 



^ + 2A = ^V A V^>(2.4) 



where is a covariant derivative, 



G^ u = R^ u - -Rg^u , Tfj, u = {p + p)u^u u - pg^ u , = <5* ■ (2.5) 

Here p(t) and p{t) are the matter density and pressure respectively, the stress-energy tensor 
corresponds to a barotropic fluid, G^ u is the Einstein tensor. 

Here it is convenient to introduce new dimensionless variables 3 : 

*(t) = Ht)/G , e(t) = d t( f>/(VA4>), (2.6) 
H(t) = H(t)/y/X = d t a/(y/Xa), p{t) = 4vrG p/A, p(t) = 4vrG p/A. (2.7) 

Here H is the Hubble parameter, and H is its dimensionless value. In these notations, 
Friedman equations for a flat universe (k = 0) in a comoving frame (u^ = [1,0,0,0]) are: 

G A = 3H 2 = ^ + l + -e 2 - 3He, (2.8) 
A <p 2 

% = 2k + 3H 2 = + 1 - -e 2 - i - 2He. (2.9) 

A q> 2 q> 

The Klein-Gordon equation (2.4) can be rewritten as: 

2p — 6p 



+ 2 = (3 + 2w) 



+ 3He 



(2.10) 



Here and below the dot denotes the derivative w.r.t. the dimensionless time i = y/Xt. 
Equations (2.8-2.10) lead to the continuity one in the form: 

+ 3H = 0, (2.11) 



P + P 

which is consistent with the equivalence principle. 

2 A here can differ from the one in AC DM theory. 

3 Here and below the present time is denoted by the subscript "0", so Go is the the current value of the 
gravitational constant. From now and further we consider the current time moment as the initial one, so 
that to = 0. Current values of cosmological parameters are taken as initial conditions. New variables lead to 
0o — l, which is convenient for further calculations. 
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3 Initial values of the model parameters 

We introduce the deceleration parameter q and the dimensionless matter density (5 for 
the initial time in the following form: 

a = - (1+q) H\ ^ ^fo-w^ aazg, (3.i) 

Combining the equations (2.8-2.10) to exclude e and <p/(j), we obtain for p = at the current 
moment to the following equation: 



w 



-i 2 9 4- 2w 

H 2 {2 -q - pz) -z\ - 2H$(3z - 1) + H*(6 - 6q - 6(3z + 4/3) = 0, z= ^^(3-2) 



This equation defines Hq as a function of f3, qo and w. In the |io| >> 1 approximation 4 
equation (3.2) yields: 



,' 2(1 + 90-/3) 

_^ ( 2-g -/3)±y ■ (3.3) 



In the GR limit (w — > oo) the second term can be neglected, so, current cosmological param- 
eters values ban be established as [33, 34]: 

Hq ~ 2.3 • 10~ 18 sec -1 , po ~ 0.27 • 10 _29 g/cm 3 (accounting for baryonic and dark matter), 
q pa —0.6. We consider a dust-filled universe, thus neglecting the pressure. We rewrite the 
above expression for the cosmological constant: 

A -+ (2 - q Q )Hl - 47rG (po - Po) ~ 11-3 • 10~ 36 sec~ 2 (3.4) 

From the lunar ranging experiment (LLR) data [35] one can extract the following limitations: 
\dtG/G\(Q) < 4 • 10 _20 sec _1 , hence |eo| is a small value: | eo ] < 0.01. When \w\ » 1 we have: 

i? ~0.68, A) -0.2, ^ ss 0.4. (3.5) 

Combining (2.8) + (2.9), multiplying the result by I/Hq anc ^ substituting (ft/cj) from (2.10), we 
obtain for the initial moment to the following expression: 



60 - 1 (2 a tf) + ttlll" 



(2-g -/3) + " o ; '- " (3.6) 

-«o 



Substituting the H value from (3.3), we obtain the order of 1/t/w 5 : 



eo J 2( -} + 0°-V (3.7) 



4 Here and below (unless otherwise noted) the arrow denotes the \w\ >> 1 approximation. 

5 When calculating the rhs of the expression (3.6) we only considered the terms of order l/y/w from (3.3), 



the last term from (3.6) was neglected due to the taken accuracy. 
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4 Dust-filled Universe solution 

First of all we consider a dust-filled Universe, i.e. p = 0. As is Ref. [1], we rewrite 
field equations using / = (/>a 3 , and take into account that the expression (2.11) leads to 
p/4> = Pofo/f- 

Considering f / f = 4>/4> + 6He + 3H + 9H 2 , we combine field equations in the following 

§[ 



way: |[(2.8)+(2.9)]+(2.10)/[6 + 4w]. This yields to: 



2 2 , r ~ r \ n 2 8 + 6^ 

/-r/ 2 (/ + Po/o) = 0, V^jj^. 
The obtained equation can be straightforwardly integrated: 

fit) 



fo 



c + E + c~/E - po , E(t) = exp(r/t) 



(4.1) 



(4.2) 



where c + and c can be easily obtained from the initial data. 
One can rewrite (2.10) as: 

2/ + 2/50/0 = (3 + 2w)(<Pa 3 ). (4.3) 
With the help of Eq. (4.2) one gets the expression for the Hubble parameter from Eq. (4.2) 6 : 



ofj = f__<t = f__ 2/p 

f </> f f(3 + 2w) 



const 



L + Arlt=L- 2(c + E-c-/E + c H ) 
fo P V f v(3 + 2w)(c+E + c~/E-p ) 

6(1 + w)(c+E- C -/E) -2c t 



i](3 + 2w)(c+E + c~/E - po) 



(4.4) 



Here c H can also be determined from initial data. 

Solving (4.2) and (4.4) for the present time, one obtains the coefficient values as: 

+ _ 1 + po e + 3H _ _ 1 + po _ e + 3H _ r/e (3 + 2w) _ e + 3H 



2 2rj 2 2i] 

resulting in the following expression for the scale factor: 

E 

,V3 



V 



- = (c+E + c-/E-po) f e*P 
a 



1 



c + E 2 + c H E -c~\ dE 
~E 



3(4 + 3w) J \c+E 2 -p E + c- 



1 



l + w 



(c+£ + c~/E - po) i+3w exp 



-2c H (A-A ) 



3(4 + 3w)VA 



(4.6) 



where 



A = 4c + c~ - pi = 1 + 2po - (3H + € ) 2 /v 



2 /„2 



H -e (l + w) 



8 + 6w L 
A(E) = arctan[(2c + £ - po)/VA} 



(4.7) 



'Note that dt = dE/(r)E). 
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Figure 1. Left hand side represents a(t) /ao according to the expression (4.6) for a dust-filled universe 
with a bounce for the following parameter values: w = —1000, q = —0.6, /3 = 0.45 for the upper 
line and (3 = 0.43653 for the lower one. Time unit is 10 10 yr. Right hand side: illustration (taken 
from NASA: http://map.gsfc.nasa.gov) for a(t)/ao in ACDM-thcory, corresponding to the Friedman 
solution. 



[c + E + c~/E- po) 4+3w exp 



2c H (A-A ) 



(4 + 3w)VA 



(4.8) 



In order to keep A positive, it is necessary to set w to be rather large (\w\ » 1) and 
negative 7 . 

In the GR case |to| — ?• oo and expression (4.6) tends to the well known Friedman solution 
with acosmological constant: 



1 E + E cr 



E r 



V3H -1 



VFr 



^3 E-E cr J cr y/3H + l 

a Fr _ (V3g + l) 2 / 3 (g-g er ) 2 / 3 
a (4£)V3 



(4.9) 
(4.10) 



It is necessary to note that in this case E = E cr , A = 0, a = 0, and the scale factor a(t) 
experiences a kink (which is absent when A > 0). The Big Bang corresponds to the moment 
ti ps -1.46A" 1 / 2 , A" 1 / 2 ss 10 10 years. 



5 Non-singular cosmology 

In BDA models, the scale factor may not vanish during it's evolution back in time unlike 
in the standard AC DM one. This is a "bounce" of the scale factor from it's minimal value 
a m 7^ 0. The bounce appears in case there is a local minimum of the scale factor greater 
than zero. The parameter phase-space for the bounce case starts from a m (E) = 0, so, it is 
possible to obtain the condition for a bounce from the equation (4.6). It has the following 
form: 



A>0. (5.1) 

7 It is also importantly to mention that the w < case in BD model opens the possibility for wormholes 
existence without energy conditions violation, see [36, 37] for details. 
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Further, the time estimation for E at the bounce is: 



E m in = y/c~/c+ . (5.2) 

The exact equality E m = w c~ /c+ is satisfied when A = 0, e.q. when the scale factor vanishes 
at the local minimum. Such a scenario excludes the initial singularity and leaves the scale 
factor regular and continuous everywhere, including during the bounce (see (4.4)). 

It is possible to estimate the numerical values of the scale factor a(t) and A at the 
bounce from the following arguments. The cosmological microwave background radiation 
(CMBR) indicates that the Universe was hot and radiation-dominated at early stages of it's 
evolution [38]. Using the expression for an adiabatic expansion ah o t/«0 = 4 • 10~ 5 [39] we can 
obtain the value of A at the bounce for \w\ >> 1 from the equation (3.7), so that: 

A « 2p a 3 m /a 3 < 2p a 3 hot /a 3 « 2.6 • 10" 14 (5.3) 

This tiny value of A can only be achieved in a nearly flat Universe, i.e. when 1 + qo — (3 ~ 0. 
This result states that the LLR bound on w remains in agreement with the cosmological one 
in BDA model, as well as with the flatness of the Universe. 

The model under consideration, as well as the usual singular cosmology 8 of [1] is not 
applicable for the hot stage of the Universe. Thus the discussed results should only be used 
to obtain initial values for the hot Universe during the evolution back in time study. 



6 Solution for a hot Universe 



The analytical study of the functions behavior near the bounce can be done only for an 
ultrarelativistic pressure. During the hot phase when p = ^p equations (2.8-2.10) lead to 
the following expression: 



H + 2H 2 



o 6 + 8u>\ 



(6.1) 



When w < —1.5, we obtain positive value of Q. 

In the FRW case (when e = and —> oo), from equation (6.1) it is possible to obtain 
expressions for the hot stage that are similar to Eqs (4.9-4.10): 



1 

71 



u + i 

U-l 



U (f) = exp 



4f 

7! 



a-hot 



1) 2 U, 



hot 



.1/4 



(u hot - iyu 



■2) 



Here a new variable f is introduced. It represents a dimensionless time, that is a measure 
from the scale factor minimum (f — t) /y/A is equal to the age of the Universe and the subscript 
"hot" corresponds to the transition from the hot stage to the cold one. 

The derivatives of the scale factor (6.2) are singular (as well as for the FRW case in 
the matter-dominated Universe — see (4.9-4.10)). Remarkably, when A — > is in the cold 
phase, the second derivative a goes to +oo at the bounce while it goes to — oo in the vicinity 
of the bounce. Therefore the Hubble function appears to be rapidly growing when A — > 0. 
The situation has to be similar to the case of the radiation-dominated Universe. Since the 
hot phase matches large values of H, it ends for a short time interval f. Therefore, when 



s The article [1] presents an analytical solution for A < 0, w > 0. 
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\w\ >> 1 and e << 1 the solution of BDA is nearly indistinguishable from the FRW one 
(except in the bounce region). 

Further, we consider the series expansion of the scale factor a(f) against f near the 
local minimum (bounce). Keeping the terms up to the fourth order, it is possible to obtain: 

a = a m + ^a m H m f 2 - ^a m 6 2 #^f 4 + .... (6.3) 

Here, H m and b are constants; H m corresponds to the second derivative of the scale factor at 
the bounce, hence it should be positive when a m > 0: H m > 0. Therefore, the equations 
for the Hubble function and it's first derivative up to second order on f are: 

H 2 = H 2 m f 2 , £= Hm{1 - b2Hmf2) -H 2 m f 2 . (6.4) 

1 + H m f 2 /2 

After substituting this into (6.1), one gets: 



l + f 2 F m (|-6 2 ) 



a a 2 _ Hm 
a ° 2 ~ (l + iLfW 



Q > (6.5) 



The last inequality is satisfied automatically (see equation (6.1)) and is valid only when 
< b 2 < 3/2. 



From (6.3), one notices that at fi = 1/y Bmb 2 the second derivative of the scale factor 
changes its sign (fi is an inflection point). So we consider an additional scale factor inflection 



point compared to the FRW case. At the time = VS/y H m b 2 = y/3fi, the first derivative 
of the scale factor changes its sign. Hence, the solutions for the hot phase and the cold one 
should be matched along the f± to T2 interval. 

When H m is large 9 and b is small, starting from the time f±, the second derivative of 
the scale factor rapidly goes to a large negative value (during the time interval of order of 
fl). Meanwhile the Hubble function remains positive (up to the moment T2). Therefore, 
along the f% to fi interval the solution for the hot phase can be matched to the cold phase 
solution. Varying the values of a m , H m and b, one could achieve a smooth connection. 

7 Conclusions 

In this paper, we have demonstrated that the Friedman solution with a cosmological 
term is a degenerated case of a more generic cosmology (for example, the BDA one as a 
ground effective approximation). 

In standard FRW cosmology, the graph of the solution of the scale factor a(t) has a 
form of a vertical line before vanishing (the bounce is possible afterwards); in the BDA case 

• with w > the graph of the scale factor a(t) vanishes with a finit slope (first derivative 
remains finite at a = 0); 



3 Here and below we compare with the unit value 
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• when w < the graph of the scale factor a(t) does not reach zero (and a bounce occurs), 
so all functions remain regular. 

Further, an adequate model with a bounce can be obtained numerically because of the 
complicated structure of the theory. Here it is importantly to note that the parameter k (we 
considered the case k = in this paper) describes the flatness type of the Universe and could 
provide a leading contribution near the bounce due to the the scale factor smallness. The 
hot phase also implies the presence of a non-vanishing pressure (we considered p = 0) which 
leads to the inability of obtaining an analytical solution for the hot phase. The presence of 
the bounce in the BDA cosmological solution allows to avoid one of the greatest problem of 
cosmology: the initial singularity. 

Finally, we would like to point out that the appearance of a bounce instead of a sin- 
gularity is ratter a common effect in gravity models when an additional scalar or tensor 
contribution is taken into account. For example, when studying the interplay between curva- 
ture and Maxwell terms in Gauss-Bonnet gravity one encounters an effect of the same nature 
when the singularity is changed by a local minimum [40]. The same effects occur in Gauss- 
Bonnet cosmology with additional fields [41]. A bounce appearance was also discovered in 
many new (sometimes exotic) models with additional terms (some examples of bounce ap- 
pearance can be found at [42-45].) and in Loop Quantum Gravity [46, 47]. So, this effect 
is rather common and natural from the mathematical point of view (changing the balance 
between different term contributions) and can be used to obtain new (stronger) estimations 
of model parameters. So, based on the condition of existence of a bounce and, so, from Eq. 
(5.3) one can put a new limit on the BD parameter w in the form 

\w\ > 10 40 , w < 0. (7.1) 

This limitation is much stronger than the existing experimental one (\w\ > 50000) and future 
developments will show the connection of these values with reality. 
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